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A comparative study of overlap and staggered fermions in QCD* 
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We perform a comparative study of the infrared properties of overlap and staggered fermions in QCD. We 
observe that the infrared spectrum of the APE/HYP improved staggered Dirac operator develops a four-fold 
near-degeneracy and is in quantitative agreement with the infrared spectrum of the overlap operator. The near- 
degeneracy allows us to identify the zero modes of the staggered operator and we find that the number of zero 
modes is in line with the topological index of the overlap operator. 



1. INTRODUCTION 

In the continuum the eigenmodes of the mass- 
less Dirac operator D are described by the eigen- 
system Dtpx = i\ilj\,X € R. Since the oper- 
ator enjoys chiral symmetry, {£',75} — 0, and 
is 75-hermitian, its spectrum is guaranteed to be 
symmetric around zero. The zero modes have a 
definite chirality, ip\^Qj5tp\=o — ±1, and are re- 
lated to the topological charge of the background 
gauge field via the Atiyah-Singer index theorem 
which provides an explicit link to the axial [/^(l)- 
anomaly and hence guarantees that the fermions 
are sensitive to topology. On the other hand, 
the low-lying non-zero modes of the Dirac opera- 
tor describe the low-energy, long-distance physics 
and are therefore referred to as the infrared (IR) 
modes. 

A naive discretisation of the fermionic Dirac 
operator describes 16 flavours of fermions instead 
of one. A 4-fold degeneracy can be removed by 
distributing the spin degrees of freedom over a 
hypercube yielding the staggered Dirac operator 
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where i]fi{x) = {—l)^"<i'^". This operator en- 
joys a remnant U{1) chiral symmetry which pro- 
tects the bare quark mass from additive renorma- 
lisation, but it still describes 4 degenerate flavours 
in the continuum. These doubler fermions are 
light and mix with each other generating a com- 
plicated particle spectrum on the lattice. More- 
over, since the operator generically has no zero 
modes, there is no (exact) index theorem and it 
is not clear, a priori, whether the operator is sen- 
sitive to topology. 

There are other conceptual issues: in order to 
generate one flavour of sea quarks in a dynami- 
cal simulation, one usually employs the so-called 
fourth-root trick, i.e., one weights configurations 
by det{D^*')^/'^ instead of det{D^^) which, how- 
ever, is hard to justify from a field theoretical 
point of view. Even if it turns out that this is 
indeed a valid procedure, the question remains 
how one should match this prescription in the va- 
lence quark sector. The crucial point in this con- 
text is therefore whether the above discretisation 
of fermions provides a consistent field theoretic 
framework or whether it is just a model of QCD, 
albeit a good one. 

One way of avoiding (at least some of) the dou- 
bling problem is to suppress the flavour changing 
interactions, i.e., to reduce the mixing of the dou- 
bler modes. This can be achieved by removing 
ultraviolet (UV) noise from the gauge fields by 
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smearing the gauge links in one way or another. 
In the foUowing we refer to this procedure as UV 
filtering or UV improvement and here we only 
note that any such procedure simply corresponds 
to a O(a^) redefinition of the staggered fermion 
action. 

Another way is to make the doublers heavy, so 
that they decouple from the theory in the contin- 
uum limit. A particular realisation of this idea is 
provided by the overlap operator jlj 

Dov = p{l + 75sign(75i:)w(-/o))) (1) 

where Z?w is the Wilson Dirac operator. The 
advantages of this operator are manifold: there 
exists an exact lattice version of the continuum 
chiral symmetry which guarantees the absence 
of an additive renormalisation of the quark mass 
and avoids mixing of operators in different chiral 
representations. Furthermore, the operator pos- 
sesses exact zero modes and allows to define an 
index theorem on the lattice implying the correct 
sensitivity of the operator to topology. 

Here we report on our results from a compar- 
ative study of the spectra of the (UV improved) 
staggered and overlap Dirac operator 0. We con- 
centrate on the IR spectra of the staggered oper- 
ator and compare it to the corresponding ones 
of the overlap operator as a benchmark. The 
question we would like to answer is whether one 
can see a 'continuum-like' behaviour in the spec- 
tra of the staggered operator similar to what has 
been observed in the Schwinger model 3 , i.e. a 
four-fold near-degeneracy and approximate zero 
modes. Then we would like to know at what lat- 
tice spacing it appears, if at all, and whether any 
deviations scale like 0(a^)-artefacts. The final 
goal, however, is to eventually test whether 

det(L',t)^/^ oc det(Dov) + 0{a^). (2) 

2. RESULTS 

We study the infrared spectra on four ensem- 
bles of quenched lattices with lattice spacings 
ranging from 0.18 to 0.07fm matched to have a 
fixed physical volume V ~ (1.12fm)'*. We use 
standard smearing techniques and we note that 
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Figure 1. Overlap and staggered eigenvalues for 
one of our quenched configurations at a ~ 0.07fm. 

for the overlap operator we use p = 1.0 through- 
out. We find that the overlap spectra are fairly 
stable under smearing for a < O.lfm. This in- 
dicates that the overlap operator separates well 
IR modes from UV fluctuations at small enough 
lattice spacings. The staggered spectra, on the 
other hand, change dramatically under repeated 
smearings. In particular, we find that without UV 
filtering there is no resemblance at all between 
the staggered and the overlap spectra. With in- 
creased UV filtering the staggered spectrum at 
large enough [3 develops an (approximate) four- 
fold degeneracy and zero modes start to sepa- 
rate out. This suggests that an (approximate) 
index theorem holds for UV filtered staggered 
fermions. Moreover we find a one-to-one corre- 
spondence between the quadruples of staggered 
eigenvalues and the ones from the overlap oper- 
ator. Occasionally, we find a mismatch of topo- 
logical charge between the index of the overlap 
and the staggered Dirac operator. In particular 
we observe that whenever the overlap operator 
has difficulties to decide what charge to assign to 
a configuration, the staggered operator will have 
problems to produce a corresponding zero mode 
- instead, the four-fold degeneracy of the would- 
be zero mode remains rather approximate. We 
assign this behaviour to the occurrence of dislo- 
cations with a size of 0(1) lattice spacings which 
sometimes are detected as topological charges and 
sometimes not. It is clear, however, that even 
the overlap operator suffers from this ambiguity 
(generically, this also happens when changing the 
parameter p in eq.JQ)) which suggests that this is 
a genuine 0{a?) effect. 

In figure n we plot the chirally rotated overlap 
eigenvalues most-left and those of the staggered 
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Figure 2. Staggered eigenvalues (geometric 
mean) vs. overlap eigenvalues at a ~ 0.07fm. 
With only 1 HYP filtering, the overlap and stag- 
gered topological charges disagree on one config- 
uration, with 7 HYP steps there is no such mis- 
match. 

operator ordered according to increasing size for 
a generic configuration with a = 0.07fm and UV 
filtering level HYPS T . After rescaling the stag- 
gered spectra with the same factor for all configu- 
rations, they show a quantitative agreement with 
the corresponding spectra of the overlap operator 
(up to mismatches in the topological charge as 
discussed above). 

In figure 121 we plot the geometric mean of the 
staggered eigenvalues against the eigenvalues of 
the corresponding overlap operator for filtering 
levels HYPl and HYP 7 on our quenched ensem- 
ble at a ~ 0.07fm. It is evident that there is 
a linear correspondence between the UV filtered 
overlap eigenvalues and the staggered eigenvalue 
quadruples at small accessible couplings and with 
only moderate UV filtering. As emphasised be- 
fore, the correspondence of the (would-be) zero 
modes suggests an index theorem for improved 
staggered fermions up to O(a^) artefacts. The 
one-to-one correspondence also implies that im- 
proved staggered fermions show the same sensi- 
tivity to topology as overlap fermions and fall into 
the same universality class of random matrix the- 
ory 0]. Indeed, this has recently been verified 
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3. CONCLUSIONS AND OUTLOOK 

We have found a four-fold degeneracy in the 
IR spectra of improved staggered Dirac opera- 
tors up to O(a^) artefacts, hence providing evi- 
dence that the rooting-procedure of the staggered 
determinant might be justified in the sense of 
eq. (O- We found a quantitative correspondence 
between the IR part of the improved staggered 
and overlap spectra in QCD on individual con- 
figurations and our data are compatible with the 
assumption that the remaining artefacts scale as 
O(a^) with moderate UV filtering and at acces- 
sible couplings. However, a number of concep- 
tual issues within the staggered fermion formu- 
lation of QCD remain. In particular it is not 
clear, whether and how a local operator D with 
det(_D) = det(Dstag)^^'' can be constructed with 
the staggered degrees of freedom. 

Finally, we note that our results could point the 
way towards unquenching the overlap Dirac oper- 
ator via a 'hybrid HMC algorithm: one could use 
an improved staggered Dirac operator along the 
molecular dynamics trajectory and accept /reject 
the final configuration at the end of the trajec- 
tory with an overlap Dirac operator. The hope is 
that in this way a correct dynamical overlap en- 
semble could be generated essentially at the cost 
of a dynamical staggered fermion simulation. 
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